Experimental evidence of conformal invariance in soap film turbulent flows 
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We present experimental evidence of statistical conformal invariance in isocontours of fluid thick- 
ness in experiments of two-dimensional turbulence using soap films. A Schlieren technique is used 
to visualize regions of the flow with constant film thickness, and association of isocontours with 
Schramm-Lowner evolution (SLE) is used to identify conformal invariance. In experiments where 
an inverse energy cascade develops, statistical evidence is consistent with such an association. The 
diffusivity of the associated one-dimensional Brownian process is close to 8/3, a value previously 
identified in isocontours of vorticity in high-resolution numerical simulations of two-dimensional 
turbulence (D. Bernard et al., Nature Phys. 2, 124, 2006). In experiments where the inverse energy 
cascade is not sufficiently developed, no statistical evidence of conformal invariance is found. 

PACS numbers: 47.27.-i, 92.60.hk, 68.15.+e, 11.25. Hf 



Two-dimensional (2D) turbulence displays important 
differences with the three-dimensional (3D) case [l|. 
While turbulent flows in 3D tend to create smaller scales 
and disorder, flows in 2D tend to self-organize and cre- 
ate long- living coherent structures 0. These differences 
are not just of academic interest, as many flows in na- 
ture are quasi-2D. This is often the result of volume 
forces that impose a preferred direction, as e.g., rota- 
tion and/or stratification in the atmosphere, and mag- 
netic fields in the interplanetary medium. In particular, 
the atmosphere can be considered to a good degree of 
approximation as a shallow layer of fluid, as most of the 
weather takes place in a thin layer (the troposphere) of 
approximately 10 km depth, while the larger horizontal 
scales are of the order of thousand kilometers [|| . 

2D turbulence is far from being completely understood. 
If injected at intermediate scales, enstrophy (the squared 
vorticity) is transferred towards smaller scales, in a so- 
called direct cascade, while energy is transferred towards 
larger scales in an inverse cascade 0, Q . Besides the di- 
rection of the cascade, and unlike its 3D counterpart, the 
energy cascade is scale invariant: while in 3D turbulence 
strong intermittent events develop that make the flow 
multi-fractal, 2D flows are absent of such intermittency 
at large scale [3, Q . 

An even stronger form of symmetry was identified in 
recent years: that of conformal invariance [(|. While 
scale-invariance is a global property, often quantified us- 
ing two-point correlation functions, conformal invariance 
is a local property involving invariance under confor- 
mal (angle-preserving) transformations. Conformal in- 
variance is used in quantum field theories, including ap- 
plications in condensed matter Q • In fluid dynamics, it is 
often used to solve Laplace equation for ideal incompress- 
ible and irrotational flows in 2D domains. That 2D tur- 
bulence (a much harder problem where non-linear effects 



are dominant) may be conformal invariant was originally 
proposed in |8(. Its confirmation requires measurement 
of multi-point correlation functions. Instead, in a recent 
work evidence of conformal invariance in numerical 
simulations of 2D turbulence was shown using Schramm- 
Lowner evolution (SLE) @. SLE applies to a special set 
of random curves without self-intersections in a 2D space, 
and it was shown that if the conformal transformation 
that maps curves in the 2D space into the real axis gives 
a one-dimensional (ID) Brownian process (BP), then the 
original curves are statistically conformal invariant. This 
result was used to show that isocontours of zero vorticity 
in simulations of 2D turbulence are conformal invariant 
(and moreover, they are generated by an SLE process) 
Q- The analysis also allowed identification of the class 
of universality to which these curves belonged, identified 
by the diffusivity re of the BP, and noted SLE K . These 
classes share common statistical properties, such as criti- 
cal exponents. More recently, the result was used to show 
that isocontours of temperature in a quasi-geostrophic 
model also were SLE, with a different re [Io| . 

It is unclear however how these results can be extended 
to more realistic configurations. Under favorable condi- 
tions, flows in nature can be considered as quasi-2D at 
best. They are not exactly 2D as they may have small 
(but finite) thickness, or translation symmetry in the out- 
of-plane direction. As a consequence, motions outside 
the plane of interest can develop. Arc conformal invari- 
ant properties still relevant in such cases? Moreover, the 
method used to detect conformal invariance in 2D tur- 
bulence relies on the vorticity, a quantity which is hard 
to measure with small uncertainties in experiments or 
observations. In many cases, vorticity is obtained af- 
ter applying a curl to particle image velocimetry (PIV) 
measurements. This results in a loss of spatial resolution 
(PIV gives, as a rule, the velocity field every six or more 
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FIG. 1: (Color online) Left: Typical raw Schlieren image of 
the soap film for configuration A. Right: Extraction of paths 
with constant thickness. The shortest paths in light gray 
(green), whose lengths are below l tI , are eliminated. 



image pixels, per virtue of the interrogation window used 
for the correlation) , and in an increase of error per virtue 
of the numerical derivatives required to compute the curl. 
In some cases, spatial resolution of the correlation tech- 
nique can be improved by using flow optic approaches 
with PIV images. However, an artificial smoothing of 
the results may still take place as a consequence of the 
regularization employed. 

In this letter, we show evidence of conformal invari- 
ance in experiments of turbulence in soap films pTH13| . 
Such experiments approximate 2D turbulence, although 
deviations arise from variations in the soap film thick- 
ness, and also because of the friction in the presence of 
rough boundari es Il4f. Indirect vorticity measurements 
have been done [l2|, LL3| an d despite these effects, its dy- 
namic is in some cases consistent with predictions for 2D 
turbulence. Also, to lowest order in the film thickness, 
thickness is advected as a passive scalar, and a strong 
correlation between thickness fluctuations and vorticity 
has been reported [3]. However, both weak and strong 
intcrmittency effects have also been observed, in the di- 
rect and the inverse cascade, attributed to finite thickness 
effects [l5| . The experiments we present here share these 
properties; a detailed description of PIV measurements 
is left for the future. 

Given these similarities and differences, it is interesting 
to know whether conformal invariance can be identified 
in the system at hand. In the experiment, grid turbu- 
lence is generated using grids with different spacing of 
wires, while the fluid flows through a vertical channel 
as a result of gravity. To avoid indirect vorticity mea- 
surements and loss of resolution associated with PIV, we 
rely for the analysis on raw images of film thickness vari- 
ation. The thickness variation in the film is visualized 
using a Schlieren technique [l6j . When employing this 
technique, gradients of film thickness act on light rays as 
density gradients do in volumetric flows. Changes in the 
thickness act as changes in the refraction index, and the 
deviation of the light rays can be used to obtain an inten- 
sity image proportional to thickness variations. The use 
of this method allows us to have images with high con- 
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FIG. 2: (Color online) (a) The zipper algorithm: the original 
path (black) is gradually conformally squeezed down onto the 
a;-axis (indicated by the gray paths and the dashed arrows). 
This defines a ID path a„(T n ) (see [13] ). (b) Cumulative 
PDF of the increments Sa T , displayed here for L = 800 pixels 
for configurations A (black dots) and B (gray or red dots); 
the dashed (blue) line corresponds to a cumulative (normal- 
ized) Gaussian distribution. The inset shows the Kolmogorov- 
Smirnov probability Pks, testing the null hypothesis that the 
increments fall into a Gaussian distribution. Black lines are 
for dataset A and long-dashed gray (red) lines for dataset B; 
the dashed (blue) line indicates a conventional threshold at 
0.1. (d) Scaling with error bars of the variance of a with 
driving time r; the best fit for configuration A (black dots) 
yields k — 2.5 ± 0.3 (solid blue line). Configuration B (gray 
or red dots) has a less clear scaling and does not pass the 
Gaussianity test. 



trast and resolution. Conformal invariance is then exam- 
ined studying isolevels of constant film thickness. As an 
illustration, we present two experimental configurations: 
one with a grid size such that a wide inverse cascade of 
energy develops, and another with a grid size such that 
only a short range of inverse cascade is available. In the 
former case, evidence of conformal invariance is found, 
while in the latter conformal invariance is not observed. 

The experimental setup is as follows. We use a z-typc 
Schlieren configuration with two mirrors of 50 cm in di- 
ameter and 3.0 m focal length. Two razor blades per- 
pendicularly disposed at the sagital and meridional focal 
planes enable us to obstruct the deviated light rays and 
to improve the contrast of the otherwise shadowgraph 
image that would form at the CCD of a monochromatic 
high-speed video camera (3000 fps with a resolution of 
512 x 512 pixels). Images are captured at a rate of 2500 
fps, with one pixel corresponding to 0.161 mm (corre- 
sponding to an area of 82.4 x 82.4 mm 2 centered in the 
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middle of the channel width). The vertical soap film 
channel has a width of 160 mm, and a length of 0.6 m 
in the region of constant width (the expansion and the 
contraction of the tunnel have the same length). In the 
first configuration (configuration A), a grid with a sepa- 
ration between wires of 5.1 mm and a width of 65 mm 
is used to create turbulence downstream. The mean flow 
velocity is 1400 mm/s, the mean film thickness is 5.6 /im, 
and the Reynolds number (based on the injection scale) 
is Re « 870. When the third order structure function is 
computed from PIV data, a range with negative energy 
flux (corresponding to an inverse cascade) can be identi- 
fied for scales between sa 30 and 110 mm. In the second 
configuration (configuration B), a grid with a separation 
between wires of 0.3 mm and a width of 45 mm is used 
instead. The mean flow velocity is 1930 mm/s, the chan- 
nel width is 160 mm, the mean film thicknes is 4.1 /im, 
and the Reynolds number is Re ~ 70. Configuration B 
has an estimated inverse cascade range (from third order 
structure functions) between w 25 and 45 mm. 

Figure [1] shows a Schlicrcn image of the soap film. In 
each analysis, 2000 snapshots of the film are used. An 
average of all snapshots is used as the background image, 
that we remove to undertake the analysis. The resulting 
image corresponds to thickness gradients, with zero value 
associated to constant thickness. For each image we then 
look at its four possible orientations, and for each one of 
those we set up an arbitrary a;-axis (see the right panel in 
Fig[T]). From each axis, we then explore all curves with no 
thickness variation until the curve goes out of the image 
or it returns to the £-axis. When the paths are not self- 
avoiding, we erase the loops similarly to how the loops 
are erased in loop erased random walks (LERWs) j9|. 
Finally, only the paths whose length are above a given 
threshold itr are kept. Here we consider In = 500 pixels. 
We thus obtain a collection of paths which are correlated. 
In order to avoid redundancies (as the flow has a typical 
time-correlation given by its turnover time), we impose 
that the paths coming from two different snapshots have 
to be separated by a certain number of snapshots nm- 
The results given below are obtained with = 8, but 
we checked that n fil = 4, 8, 12, 16, and 20 give the same 
results. For ngi > 20, similar results are obtained al- 
though with worse statistical accuracy as less curves are 
preserved (nsi = 20 roughly corresponds to the turnover 
time of the eddies at the injection scale in configuration 
A). With nfii = 8, the number of paths above 500 pixels 
is 1934 for configuration A, and 2216 for configuration B. 

These curves are projected with a conformal map into 
ID paths (the driving functions, which we will label 
a(r)), as done in |(|. The zipper algorithm is used to 
obtain the ID paths; see and Fig. [21(a). The para- 
metric variable r (the "time" ) here only labels the succes- 
sive points of the ID paths, and should not be confused 
with the actual time corresponding to the different snap- 
shots in the experiment. For increasing r, all curves in 
the 2D plane are generated in the direction that keeps 
positive vorticity to the right. Following Schramm's re- 
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FIG. 3: (Color online) Probability that a path keeps a given 
point za to its right. The dashed (blue) line indicates the 
SLE prediction given by Eq. JTJ with k = 8/3. Black dots are 
for dataset A and gray (red) dots for dataset B. In the latter 
a larger departure is observed from the predicted behavior. 
The inset shows a sketch of a 2D path that keeps the given 
point zo to its left, and indicates how the angle 9 — arg(zo) 
is defined. 

suits [§] , in order to check whether the measure of these 
curves is conformal invariant, we have to test whether 
these driving functions are likely to be seen as Brown- 
ian ID paths. In other words, the paths must be of the 
form a(r) = y / K_B(r), where n is a diffusion coefficient, 
and B(t) is a standard Brownian motion. To test this, 
we first perform a Kolmogorov-Smirnov test, testing the 
null hypothesis that the increments of the driving func- 
tions fall into a Gaussian distribution (whose mean and 
variance are data driven). To compute the increments, a 
length L is picked, and all driving functions with length 
larger than L are selected. Then, the minimum driv- 
ing time t for which L is reached in all these functions 
is found, and the increments are defined as 6a T = a(r) 
with r the minimum driving time for the given L. As 
an example, the cumulative probability density function 
(PDF) of the increments for datasets A and B and for 
L = 800 pixels is shown in Fig. [2Kb). The pks values 
of the Kolmogorov-Smirnov test for different minimum 
driving times r (i.e., different lengths L) are shown in 
the inset of Fig. [2Kb) . For configuration A the null hy- 
pothesis is rejected. Finally, we study the dependence of 
the variance of a as a function of the driving time r for 
datasets A and B to evaluate k; see Fig. [21c). For con- 
figuration A we obtain k ss 2.6; configuration B presents 
slightly larger deviations from the linear behavior and 
does not pass the Gaussianity test. 

The procedure described above was also tested on syn- 
thetic 2D paths generated from LERWs, and from self 
avoiding walks (SAWs). The former are generated with 
a 2D Brownian motion in the upper plane whose loops are 
successively erased, while the latter are generated with 
a Markov chain Monte Carlo algorithm (the "pivot algo- 
rithm" [HI)- in each case, 1000 paths arc used. Both 
synthetic datasets pass the Gaussianity test and have 
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variance re = 2.0 ± 0.2 (LERWs) and 2.9 ± 0.3 (SAWs). 

The driving functions for configuration A show good 
Pks values as well as good scaling of the variance with 
t. The analysis thus indicates that the isolines used are 
good SLE candidates. The value found is compatible 
with k = 8/3, reported previously in numerical simu- 
lations @, under the hypothesis that we are measuring 
here the properties of the envelope of the paths, thus 
accessing to the dual. This value corresponds to the con- 
tinuum limit of SAWs. Configuration B, with a smaller 
Reynolds number and a narrower inverse energy cascade 
range, may not be sufficiently turbulent for conformal 
invariance to develop and be identified. 

To confirm these findings, we now consider two prop- 
erties that isocontours should exhibit if they are confor- 
mal invariant. First, we calculate the fractal dimension 
df of each of the paths, and study the distribution of 
k = 1 + df /8. For the two synthetic datasets mentioned 
before, the distributions of k based on the fractality of the 
2D paths have mean 1.9 and variance 0.3 in the LERWs, 
and 2.4 ± 0.5 in the SAWs. For datasct A, k = 2.4 ± 0.4, 
while for datasct B, k = 2.3 ± 0.6. In other words, frac- 
tality (as given by a measure of df) obtains in both cases 
A and B, and we observe good agreement between k and 
the value derived from the distribution of the fractal di- 
mensions of the 2D paths. 

Second, we verify a non-trivial result predicted by the 
SLE theory. If the curves are conformal invariant, then 
the probability that a path keeps a given point zq to its 
right depends only on k and 9 = arg(zo), namely fl7| 
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where L is the Gamma function, and 2*1 is the Gauss 
hypcrgcomctric function. This is a good example of non- 
trivial exact predictions that can be obtained if 2D or 
quasi-2D turbulence is conformal invariant (see for 
other examples). Results are shown in Fig. [3j Although 
both configurations show good agreement with the pre- 
diction, the agreement is better for dataset A. 



The results indicate that the soap film flows studied 
here exhibit a behavior consistent with conformal invari- 
ance when an inverse energy cascade is sufficiently de- 
veloped, at least to the point that SLE predictions are 
in good agreement with the experimental results. Note 
that an experimental result for SLE in a different context, 
that of the rugosity of tungsten oxide surfaces, has been 
reported recently [19|, identified in their case with the 
Ising model (k = 3). Given the role played by 2D and 
quasi-2D flows as simple (although incomplete) models 
of atmospheric an oceanic turbulence, it is to be hoped 
that studies which bring together fluid dynamics and crit- 
ical phenomena, together with the link to conformal field 
theory, will shed new light on the statistical properties of 
such flows. 
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